Perimeter Length and Form Factor of Two-Dimensional Polymer Melts 
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Self-avoiding polymers in two-dimensional (d = 2) melts are known to adopt compact configu- 
rations of typical size R(N) ~ N 1 ' d with N being the chain length. Using molecular dynamics 
simulations we show that the irregular shapes of these chains are characterized by a perimeter 
length L(N) ~ R(N) dp of fractal dimension d p = d - 6 2 = 5/4 with 9 2 = 3/4 being a well- 
known contact exponent. Due to the self-similar structure of the chains, compactness and perimeter 
fractality repeat for subchains of all arc-lengths s down to a few monomers. The Kratky repre- 
sentation of the intramolecular form factor F(q) reveals a strong non-monotonous behavior with 
q 2 F(q) ~ l/(qN 1/d ) @2 in the intermediate regime of the wavevector q. Measuring the scattering of 
labeled subchains the form factor may allow to test our predictions in real experiments. 

PACS numbers: 61.25.hk 
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Introduction. It is well known that linear polymers 
in two dimensions (2d) adopt compact and segregated 
conformations at high densities [H, 13, [3, 0, H, Q ■ This 
is expected to apply not only on the scale of the total 
chain of N monomers but also to subchains comprising 
s monomers, at least as long as the segments are not 
too small (1 <C s < N). The typical size R(s) of a chain 
segment should thus scale as R(s) ~ s v with an exponent 
v = 1/d set by the spatial dimension d = 2. Compactness 
does obviously not imply Gaussian chain statistics 0, Q 
nor does segregation of chains and chain segments impose 
disk-like shapes minimizing the average perimeter length 
L(s) of chain segments. The boundaries of chains and of 
chain segments are in fact found to be highly irregular 
as revealed by the snapshot presented in Fig. [TJ Using 
scaling arguments and molecular dynamics simulations 
we show below that these perimeters are fractal, scaling 
as 
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<, p = d — 02 = 5/4 > 1 being the fractal line dimen- 
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sion. Our work is based on the pioneering work by Du- 
plantier who predicted a contact exponent 62 = 3/4 [2j 
characterizing the intrachain segmental size distribution. 
In contrast to many other possibilities to characterize nu- 
merically the compact chain conformations the perimeter 
length can be related to the intrachain form factor F(q) 
making it accessible experimentally, at least in principle, 
by means of small-angle scattering experiments 0, HI • 

We recall first the computational model used for this 
study and confirm then numerically the scaling of R(s) 
and L(s) suggested above. The analysis of intrachain 
properties such as the segmental size distribution, the 
bond-bond correlations, and the intrachain form factor 
will allow us to demonstrate Eq. (fT]). We conclude by 
discussing consequences for the dynamics of 2d melts. 

Computational issues. Our numerical results are 
obtained by molecular dynamics (MD) simulations of 
monodisperse linear chains at high densities. Our coarse- 




FIG. 1: (color online). A snapshot of a melt configuration 
with chain length = 1024 and monomer density p = 7/8 
can be seen in the left panel. Only the perimeter monomers 
interacting with other chains are indicated. The numbers 
refer to an arbitrary chain index used for computational pur- 
poses. The chains are compact, i.e., they fill space densely, but 
compactness does not imply a disk-like shape (see, e.g., chains 
3 and 8). The main figure presents the end-to-end distance 
R{s) = (r 2 ) 1 / 2 , the radius of gyration R s (s), and the mean 
perimeter length L(s) of segments containing s = m — n + 1 
monomers as indicated by the sketch on the right. The full 
symbols refer to overall chain properties (s = N). The solid 
lines confirm the exponent v = 1/2 for the segment size, the 
dashed line confirms the scaling of L(s) suggested by Eq. {T}. 



grained polymer model Hamiltonian is essentially iden- 
tical to the well established Kremer-Grest (KG) bead- 
spring model 0, [Tot 11 1 where the excluded volume inter- 
action among monomers is mimicked by a purely repul- 
sive Lennard- Jones potential and the chain connectivity 
is assured by harmonic springs calibrated to the "finite 
extendible nonlinear elastic" (FENE) springs of the KG 
model. We focus in this work on melts of density p = 7/8 
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at temperature ksT = 1 with chain lengths ranging up 
to N = 2048 using a periodic simulation box of linear 
length 335.2 containing 98304 monomers 

Mean segment size and perimeter length. The main 
part of Fig. Q] presents the typical size and perimeter 
length of a chain segment between the monomers n and 
m = n + s — 1 as indicated by the sketch. Following 
[l3L [I4I ] we average over all pairs (n, m) possible in a 
chain of length TV. Averaging only over segments at the 
curvilinear chain center slightly reduces chain end effects, 
however the difference is negligible for the larger chains, 
TV > 256, we will focus on. Open symbols refer to seg- 
ments of length s < TV of chains of length TV = 1024, full 
symbols to total chain properties (s = TV). The segment 
size may by characterized by either the second moment 
R 2 {s) — (r 2 ) of the end-to-end vector r of the segment 
(squares) or by its radius of gyration i? 2 (s) (triangles) 
[ia ]. In agreement with various numerical studies 
the presented data confirms that the chains are compact, 
i.e. v = 1/2 (solid lines), on all scales s [l6j . A perime- 
ter monomer of a chain segment is defined as a monomer 
being within a distance 1.2 to a monomer not belonging 
to the same chain segment 17J. The mean number L(s) 
of these perimeter monomers increases with a power-law 
exponent 1 — 1/Q2 = 5/8 (dashed line) which is in perfect 
agreement with Eq. ([1]) , and this holds again on all scales 
for arbitrary segment lengths provided that the segment 
is sufficiently large (s > 50). We demonstrate in the 
following where the suggested scaling stems from. 

Segment size distributions and contact exponents. Ob- 
viously, the mere fact that the exponent v is the same 
in 2d and 3d does not imply that 2d melts are Gaussian 
0. This can be directly seen from the different proba- 
bility distributions of chain segment vectors r = r m — v n 
presented in Fig. [2l To simplify the plot we focus on the 
two longest chains simulated, TV = 1024 and TV = 2048. 
Go (r, TV) characterizes the distribution of the total chain 
end-to-end vector (n = 1, m = TV), Gi(r, TV) the distance 
between a chain end and the monomer in the middle of 
the chain (n = 1, m = TV/2) and G2(r, TV) the distribu- 
tion of an inner segment vector between the monomers 
n = N/4 and m = 3TV/4. In addition, we indicate the 
segmental size distribution G(r, s) averaging over all pairs 
(n, m) which has been used recently to characterize de- 
viations from ideal chain behavior in 3d melts 14j . All 
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data for different N and s collapse on three distinct mas- 
ter curves if the axes are made dimensionless using the 
second moment R 2 of the respective distribution as in- 
dicated in the figure. The only relevant length scale is 
thus the typical size of the segment itself. The distri- 
butions are not monotonous and are thus qualitatively 
different from the Gaussian (thin lines) expected for ran- 
dom walks. In agreement with Duplantier Q we find 
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FIG. 2: (color online). Scaling plot of various distributions 
of the end-to-end vector r = r m — r n of chain segments of 
chains of length N = 1024 (stars) and TV = 2048 (triangle): 
G (r,N) for n = 1 and m = TV, Gi(r,N) for n = 1 and 
m = TV/2, G 2 (r,TV) for n = N/4 and m = 3TV/4. G (r,TV) 
and Gi (r, TV) are shifted vertically for clarity. The segmental 
size distribution G(r, s) averaging over all pairs (n, m) given 
for TV = 1024 with s = 256 (squares) and s = 512 (spheres) 
scales as G<z(r, TV). The thin lines indicate the Gaussian distri- 
bution y — exp(—x 2 ) /it expected for ideal chains in 2d. The 
power laws y « x @i (dashed lines) observed for 1 < 1 have 
been predicted by Duplantier Q]. The Redner-des Cloizeaux 
formula for G2 (r, TV) is indicated by a solid line. 



exponents Oo = 3/8, ©1 = 1/2 and O2 = 3/4 (dashed 
lines) describing the small- a: limit where the universal 
functions fi(x) become constant. Especially the largest 
of these exponents, O2, is clearly visible. The contact 
probability for two monomers of a chain in a 2d melt is 
thus strongly suppressed compared to ideal chain statis- 
tics (©0 = 61 = 62 = 0). As can be seen, the rescaled 
distributions G(r,s) and G2(r, TV) become identical for 
intermediate segment length, 1 C s C TV. (Obviously, 
G(r, s) « Go(r, TV) for very large segments s — > TV.) It is 
for this reason that the exponent ©2 is the most impor- 
tant one for asymptotically long chains where chain end 
effects can be neglected. The rescaled distributions show 
exponential cut-offs for large distances. The Redner-des 
Cloizeaux formula fl8j is a useful interpolation formula 
which supposes that fi(x) — Ci exp(— KiX 2 ) with con- 
stants Ci and Ki = 1 + ©i/2 imposed by the normal- 
ization and the second moment of the distributions [l^ | . 
This formula is by no means rigorous but yields reason- 
able parameter free fits as shown for f^ix). 

Bond-bond correlations. The bond-bond correlation 
function -Pi(s) = (e„ • e m ) (ej denoting the normalized 
bond vector connecting the monomers i and i + 1) has 
been shown to be of particular interest for characterizing 
the deviations from random walk statistics in 3d polymer 
14j . The reason for this is that Pi(s) is propor- 
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with x — r/Ri being the scaling variable and the contact tional to the second derivative of the segment size R(sY 
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FIG. 3: (color online). Bond-bond correlation functions 
Pi(s) = (e„ ■ e m ) and P2(s) = ((e„ ■ e m ) 2 ) — 1/2 vs. curvilin- 
ear distance s — 1 between normalized bond vectors e n and 
e m . The first Legendre polynomial Pi(s) (inset) shows an 
anti-correlation at s a 12. The second Legendre polynomial 
decays over two orders in magnitude as a power law 
(dashed line) with an exponent 1 + vQ-z — 11/8 in agreement 
with the return probability calculated from Eq. @. 



with respect to segment length s so that small deviations 
from the asymptotic exponent 1v = 1 are emphasized 
14j . As can be seen from the inset in Fig. [3] deviations 
of this kind are small for large s and may be neglected 
for the present study. The main effect visible is an anti- 
correlation at s a 12 due to the backfolding of the chain 
contour which can be directly seen from the snapshot. 
Conceptually more important is the fact that the second 
Legendre polynomial P2(s) = ((e„ • e m ) 2 ) — 1/2 reveals 
a clear power law behavior over two orders of magnitude 
in s (dashed line). The power law is due to the align- 
ment of two bonds if they are sufficiently close, i.e. the 
exponent measures the return probability after s steps. 
It follows from Eq. ([2]) that for 1 <C s <C N this is given 



by linv^o G(r, s) ~ 1/s 
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The agreement 



of the data with this exponent is excellent and provides 
an independent confirmation of 02 =3/4. 

Intrachain form factor. Neither segmental size dis- 
tribution nor bond-bond correlation function are read- 
ily accessible experimentally. It is thus important that 
02 should be measurable — at least in principle - 
from an analysis of the intrachain form factor F(q) = 

W J2n,m=i( ex P ' ( r ™ ~ r ™)D- Tne reason for this is 
that the form factor can be expressed by the Fourier 
transform G(q, s) = (cxp [iq ■ (r„ - r m )]) nm of the seg- 
mental size distribution G(r, s) with q being the wavevec- 
tor: F(q) = jjj ds(N - s)G(q, s). Assuming G(r, s) a 
Gi (r, TV) and using the Redner-des Cloizeaux approxima- 
tion (Eq. ([2|)) this yields a lengthy analytic formula (no t 
given) which is represented by the solid line in Fig. [3] 20]. 
For wavevectors corresponding to the power-law regime 
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x = qR g (N) 



FIG. 4: (color online). Kratky representation of the in- 
tramolecular form factor F(q) as a function of x = qR g (N) 
for different chain length N using the same symbols as in 
Fig. [3] The Debye formula (thin line) corresponds to a chain 
length independent plateau for i > 1. By contrast to this, 
a strong non-monotonous behavior is revealed by our data 
which approaches with increasing N a power law exponent 
— 02 = —3/4 (dashed line) corresponding to a compact object 
of fractal line dimension d p = d — 02 = 5/4. Also included is 
the Porod scattering expected for a compact 2d object with 
smooth perimeter (dash-dotted line) and the Fourier trans- 
form of the Redner-des Cloizeaux approximation (solid line). 
The increase of the scattering for large q (Bragg peak) is due 
to the packing of the beads on local scale. Only in this limit 
does F(q) become chain length independent. 



of Eq. ((2]) this reduces to the simple power law 

F(q)/N a 1.98/ \qR e (N)) d+e * (3) 

indicated by the dashed line. Note that the above scal- 
ing is a direct consequence of G(r, s) ~ r 02 and does 
not rely on the Redner-des Cloizeaux approximation. 
We rescale the wavevector with the measured radius 
of gyration R g (N) (presented in Fig. Q]) to collapse all 
data in the Guinier regime for small x = qR g (N) <C 1 
and use a Kratky representation for the vertical axis 
y = (F(q)/N)x 2 . While y becomes constant and inde- 
pendent of chain length for x 3> 1 for Gaussian chains (as 
shown by the Debye formula indicated) [l|, we observe 
over a decade in i a striking non-monotonous behavior. 
Our data suggests that Eq. ([3]) is approached systemati- 
cally with increasing chain length — the central numeri- 
cal result presented in this paper. 

Identification of Q2 and the fractal line dimension. 
The preceding discussion focused exclusively on intra- 
chain properties. Since 2d chains are compact (Fig. [1]) 
only monomers on the chain perimeter interacting with 
monomers from other chains can contribute to the scat- 
tering. Quite generally, the scattering intensity NF{q) of 
compact objects becomes proportional to the mean "sur- 
face" L(N) ~ R(N) d e for qR(N) <C 1 which implies the 
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generalized Porod law @, ED, [H| 

NF(q) « N 2 / (qR(N)) 2d - dp . (4) 

For a 2d object with smooth perimeter (d p = 1) this cor- 
responds to the classical Porod scattering F(q) ~ 1/q 3 
represented by the dash-dotted line in Fig. |U Compar- 
ing Eq. (|U) with Eq. shows that 2d melts are char- 
acterized by a fractal line dimension d p = d — 02 and 
demonstrates finally the scaling of the perimeter length 
L(N) postulated in the Introduction and verified nu- 
merically in Fig. [U By labeling only the monomers of 
sub-chains (which corresponds to a scattering amplitude 
sF(q) ~ L(s) ~ R(s) dp ) the above argument is read- 
ily generalized to the perimeter length L(s) of arbitrary 
segment length s < N |23j |. 

Summary. Investigating various static properties of 
linear polymer melts in 2d we have demonstrated that 
the compact chains and chain segments (y = 1/2) are 
characterized by a fractal perimeter L(s) of line dimen- 
sion dp = d—®2 = 5/4. As may be seen from Fig. [4] com- 
putationally very demanding systems with chain length 
N > 10 3 are required, and have thus been simulated, 
to put to the test the suggested scaling behavior. Our 
results may be verified experimentally from the scaling 
of the intrachain form factor F(q) whose Kratky repre- 
sentation is predicted to reveal a strong non-monotonous 
behavior. This should also hold in semidilute solutions 
provided that the chains are long enough '12J. Interest- 
ingly, the fractality of the perimeter precludes a finite 
line tension. Thus, shape fluctuations of the segments 
are not suppressed exponentially [H, but may occur in 
an "amoeba-like" fashion by advancing and retracting 
"lobes". According to a recent suggestion Q the relax- 
ation time t(s) of a chain segment may, hence, scale as 
t(s) ~ i(s) 3 ~ s 15 / 8 rather than as s 2 as in the stan- 
dard Rouse model [15| . Clearly, as Gaussian chain statis- 
tics is inappropriate to describe conformational proper- 
ties of 2d melts, there is no reason why a model based on 
this statistics should allow to describe, e.g., the perime- 
ter length fluctuations. Since the latter property is in 
principle accessible experimentally from the dynamical 
intrachain structure factor F(q, t) 0, EH this is an im- 
portant issue we are currently investigating. 
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